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Abstract 

This paper describes preliminary results on the inclusive process b ^ obtained fi'om 20.7 fb^^ 
of data recorded with the BABAR detector during 1999-2000. Signal event yields are found from a 
combination of twelve exclusive decay channels after subtracting continuum and BB backgrounds. 
Cross-feed from incorrectly reconstructed 6 — > S7 events is also removed. Branching fractions in 
bins of hadronic mass are calculated using corrected Monte Carlo signal efficiencies; this is equivalent 
to measuring the gamma energy spectrum. We measure the first moment of the gamma energy spec- 
trum constraining the HQET parameter A = 0.37 ±0.09 (stat) ±0.07 (syst) ±0.W (model) GeV/c^. 
A fit to the hadronic mass spectrum gives B{b — > 57) = (4.3ib0.5 (stat) ±0.8 (syst) ±1.3 (model)) ■ 
10~^ for the inclusive branching fraction. We also constrain the HQET parameter Ai. 
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1 Introduction 



The b ^ S'j transition proceeds by an effective flavor-changing neutral current. The Standard 
Model prediction, in which the t penguin loop gives the largest contribution, has been calculated 
to next -to-leading order [Q] . Measurements of the inclusive branching fraction have been used to 
constrain new physics contributions to the decay amplitude 0. 

Since 6 ^ 57 is a two-body decay process, the photon energy, E-y, in the B rest frame, is related 
to the recoil hadronic mass, Mffad: by: 



2Mb 



In our analysis we fully reconstruct twelve exclusive h ^ decays, and use the hadronic mass 
spectrum to measure the photon energy spectrum in the B rest frame. Note that in a fully inclusive 
analysis [Q, where only the photon is measured, the gamma energy is smeared by the resolution 
of the calorimeter and by the motion of the B meson in the T(45') rest frame. With our semi- 
exclusive approach, the gamma energy resolution depends only on the Muad resolution, which is a 
few MeV/c^. 

We fit the Muad spectrum with the model proposed by Kagan and Neubert Q, which predicts 
the shape of the gamma energy spectrum using heavy quark effective theory (HQET). There are 
two main parameters, an effective h quark mass, m^, and Ai, which is related to the kinetic energy 
of the b quark in the B meson. The model introduces the known K*^ contribution using local 
quark-hadron duality to convert the portion of the spectrum below a cutoff mass into a i('*(892) 
Breit-Wigner shape. The cutoff mass is a free parameter in the model, but is expected to be 
between 1.0 and 1.2 GeV/c^. Above this cutoff the inclusive model does not explicitly include any 
higher resonances. 

We find a correlated band of allowed values for the parameters mf, and Ai, and a cutoff mass of 
about 1.11 GeV/c^. By measuring the first moment of the photon energy spectrum we constrain 
the HQET parameter A and hence m;,. A fit to the Muad spectrum with this constraint applied 
to rrib gives the inclusive branching fraction and a range of values for Ai. There is currently much 
interest in the HQET parameters, since they are needed to extract Vub and Vcb from semileptonic 
B decays @. 



2 The BaBAR detector and dataset 

This paper describes preliminary results obtained from 20.7 fb^^ of data recorded with the BABAR 
detector at the PEP-II asymmetric e"'"e~ storage ring during 1999-2000. The data correspond to a 
total of (22.74 ± 0.36) • 10*^ BB pairs collected on the T{AS) resonance. 

The BABAR detector is described in detail elsewhere |^. We briefly summarize the detector 
systems most relevant to the current paper. The BABAR detector contains a five-layer silicon vertex 
tracker (SVT) and a forty-layer drift chamber (DCH) situated in a 1.5 T solenoidal magnetic field. 
These devices detect charged particles and measure their momentum and ionization energy loss 
(dE/dx). The transverse momentum resolution is CTpt/pt = (0.13±0.01)%-p( + (0.45±0.03)%, with 
Pt in GeV/c. Photons are detected in a CsI(Tl) crystal electromagnetic calorimeter (EMC). The 
EMC detects photons with energies as low as 20 MeV. The nominal EMC resolution for photons 
and electrons is oe/E = 2.3%/i?^/^ 1.9%, with E in GeV. The charged particle identification 



8 



(PID) combines dE/dx measurements in the SVT and DCH with particle velocity measurements, 
obtained with an internally reflecting ring-imaging Cherenkov detector (DIRC) of quartz bars 
surrounding the DCH. 

The interactions of particles traversing the detector are simulated using GEANT4 Q . Variations 
in detector conditions and beam-induced backgrounds are taken into account in these simulations. 
Signal and generic background samples are used to study the effect of the event selection criteria 
and to estimate the backgrounds. The generic background simulations consist of e^e~ qq 
{q = u, d, s, c) and BB events, where the b —>■ signal events have been removed. 

3 Event Selection 

We reconstruct 12 final state^ formed from a high-energy photon, a charged or neutral kaon, and 
1-3 pions, one of which can be neutral: 

These make up about 50% of the total b ^ rate for Mjjad < 2.4 GeV/c^. 

We select events with at least one photon with an energy 1.5 < E* < 3.5 GeV in the e~^e~ 
centre-of-mass (CM) frame. If the photon can be combined with another photon with energy 

> 50 MeV (250 MeV) to make a 7r° (r/) with a mass between 115 MeV/c^ (508 MeV/c^) and 
155MeV/c2 (588 MeV/c^), it is vetoed. 

In the hadronic system, a tt^ is made from two photons, with a minimum photon energy of 
30 MeV, and a minimum energy of 200 MeV. Well reconstructed tracks are used for charged 
pions, but charged kaons must satisfy additional kaon identification criteria. Candidates for Ks — >■ 
TT+vr" decays are selected by requiring a mass between 489 and 507 MeV/c^, and a decay length 

> 2 mm. These kaon selections reduce any b ^ dj background contributions to a negligible level. 

The mass of the hadronic system is restricted to: 

0.6 < Mnad < 2.4 GeV/c^ 

where the upper limit corresponds to a requirement on E^ > 2.094 GeV. Most of the background 
is accounted for by events in which the photon comes from initial state radiation (ISR) or from 
TT^/r) decays that survive the veto. We suppress these by requiring cos 9t* < 0.7, where cos^r* is 
the angle in the CM system between the high-energy photon candidate and the thrust axis of the 
rest of the event calculated after removing the B candidate. Further suppression is achieved by 
combining several event shape variables in a Fisher discriminant. The variables include the ratio 
of the second and zeroth Fox Wolfram moments computed in the CM and in the photon recoil 
system, the CM direction of the B candidate, and energy flow in a 20° cone along the photon 
direction and in a 40° cone opposite the photon. 

The hadronic system is combined with a high-energy photon, and accepted as a i? candidate if 
the energy in the CM system, E'^, is within 150 MeV of the CM beam energy, E^^^^^, : 

\AE\ = \E%-E;,,.J<150 MeV (2) 
^Charge conjugate modes are assumed throughout this paper. 
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and the beam-energy substituted mass rriES- 



"^ES = 7^^^^ > 5.21 GeV/c' (3) 

where is the B momentum in the CM frame. 

If there is more than one B candidate left in an event, we choose the candidate with the 
minimum | Ai?]. Then, we rescale the measured photon energy to give AE = GeV, and recalculate 
rriES- This corrects for energy leakage and improves the mss resolution, particularly for candidates 
without a TT^. 



4 Background Subtraction 

The background is mainly from continuum, although BB backgrounds become significant at high 
^Had values. The continuum background mss distribution can be fit by an ARGUS |1C] function. 
When divided into bins in M^adi the generic continuum Monte Carlo samples show a slight depen- 
dence of the shape parameter of the ARGUS function on Mnad, which we take into account in our 
fits. 

The BB background is mainly composed of final states with high multiplicity and high hadronic 
mass, where the photon is the daughter of a vr*^ or that decayed asymmetrically. We fit the generic 
BB Monte Carlo samples with an ARGUS function for the non-peaking background and a Crystal 
Ball function [11| for the peaking background, such as i? ^ D^*^p~ where only a low-energy 



photon is missing. The shape parameters of the ARGUS functions for the continuum and non- 
peaking BB backgrounds are sufficiently similar that it is possible to use a combined ARGUS 
function to describe the sum of the two backgrounds. 

Figure || shows the continuum and non-peaking BB backgrounds for mss > 5.21 GeV/c^, in 
bins of hadronic mass. The level of BB background is only significant for Muad > 1.8 GeV/c^. 
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Figure 1: Expected continuum and BB background events in 20.7 fb ^ from Monte Carlo (MC) 
simulation as a function of the hadronic mass. 



We treat b ^ s'j events as "cross-feed" background when they are reconstructed with the correct 
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photon, but with the wrong hadronic system, since this gives an incorrect measurement of Mffad- 
Cross-feed events come from three sources, which are, in decreasing order of importance: 

• Final states that are not considered in the analysis, e.g. K^, > 2Tr^, > 5-body hadronic final 
states, and events with Mnad > 2.4 GeV/c^; 

• Final states which are considered, but in which one of the final state particles is not detected; 
and 

• Multiple candidate events where the wrong candidate is chosen. 

The cross-feed tit-es distribution is fit with an ARGUS function plus the Crystal Ball function to 
allow for a possible peaking contribution. The cross-feed ARGUS function rises towards the signal 
region. Note that the level of cross-feed background is proportional to the inclusive branching 
fraction. The peaking component of the cross-feed background contributes just a few events. 

Fits to the ttt-es distributions of data events are performed in each of the nine Muad bins 
using a Crystal Ball function for the contributions that peak in the signal region and two ARGUS 
functions for the non-peaking backgrounds. The plots of the fits to the eight Muad bins from 0.8 
to 2.4 GeV/c^ can be seen in Figure ^. 

Table |^ lists the peaking yields for the sum of all final states in data, BB Monte Carlo and 
cross-feed, normalized to our measured inclusive branching fraction. 
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data 


BB 


cross-feed 


0.6-0.8 


10 




6 








0.0 


± 0.1 


0.8-1.0 


120 


± 


13 








4.0 


± 0.5 


1.0-1.2 


27 




8 








2.8 


± 1.0 


1.2-1.4 


70 




15 








9.2 


± 3.2 


1.4-1.6 


77 




14 








7.2 


± 4.4 


1.6-1.8 


49 




16 








5.0 


± 2.9 


1.8-2.0 


53 




14 


12 


± 


7 


7.3 


± 4.1 


2.0-2.2 


39 




17 


16 




12 


-2.0 


± 6.2 


2.2-2.4 


41 




16 


30 




15 


6.7 


± 2.2 



Table 1: The data, BB and cross-feed yields for all final states in the nine Muad bins with their 
statistical errors. 



5 Signal Efficiency 

The signal efficiency is determined from Monte Carlo samples using the yield for true signal events 
after cross-feed has been removed. In the first two Mfjad bins, 0.6-1.0 GeV/c^, the signal is modeled 
by the four exclusive K*^ modes. In the range 1.0-2.4 GeV/c^, we use the inclusive Xs"f model of 



Kagan and Neubert with rrib = 4.65 GeV/c^. This inclusive model uses non-resonant JETSET [12| 
fragmentation of the Xs to produce the hadronic final states. Changing the modeling of the 
1.0-1.2 GeV/c^ bin to -fC*7, and varying the parameters of the inclusive model form part of our 
systematic studies. 
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Figure 2: The fits to the data in bins of Muad between 0.8 and 2.4 GeV/c^. 
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The overall efficiency in each Mfjad bin is an average over the individual final states, weighted by 
the fractional contribution of that final state to the total in that Mfjad bin. These efficiencies have 
to be corrected by ~ 10% for the differences between data and Monte Carlo detection efficiencies for 
7r°, 7 and Ks- There are significant differences between the efficiencies for the individual 
final states and a strong dependence of the overall efficiency on Mnad- To illustrate this, the 
efficiencies can be understood in terms of a l/Mnad dependence of the efficiency for each final 
state, a 60% efficiency ratio between vr'^ and 7r+ final states, a 55% efficiency ratio between Ks and 
final states, and a 60% efficiency factor between 3 and 2-body, and between 4 and 3-body 
hadronic systems. 

To check the fractions of each final state in the non-resonant Monte Carlo, we divide the 
sample above Mfjad = 1-0 GeV/c^ into three sets of final state categories: with and without a vr*^, 
with and without a Kg^ and 2, 3, or 4-body hadronic systems. A check of the KgjK^ samples 
shows that the data and Monte Carlo values agree, and the ratio Ks/K~^ is consistent with 0.5 
as expected by isospin arguments. However, we find that the fraction of vr*^ final states needs to 
be increased by a factor 1.5 and the fraction of 2-body final states decreased by a factor 0.4 in 
order to obtain a reasonable agreement between non-resonant Monte Carlo simulation and data. 
When these adjustments are applied to the fractions of final states generated by JETSET, they 
lead to a correction of ~ 15% in the overall efficiency. Table |2| presents the signal efficiencies after 
all corrections have been applied. 



MHad [GeV/c^ ] 


Efficiency 


0.6-0.8 


0.133 




0.009 


0.8-1.0 


0.094 




0.002 


1.0-1.2 


0.058 




0.004 


1.2-1.4 


0.045 




0.003 


1.4-1.6 


0.040 




0.002 


1.6-1.8 


0.034 




0.002 


1.8-2.0 


0.022 




0.002 


2.0-2.2 


0.022 




0.003 


2.2-2.4 


0.015 




0.003 



Table 2: The corrected signal efficiencies in bins of M^jad with their statistical errors. 

6 Branching Fractions 

Partial branching fractions (AS) are calculated for each Mfjad bin for the sum of the twelve final 
states reconstructed in this analysis. The calculation uses the data and BB yields from Table |l| 
and the signal efficiencies from Table The peaking cross-feed background has to be normalized 
to our measured inclusive branching fraction (B). We take care of this dependence of the peaking 
cross-feed background by including effective terms that modify the signal efficiencies in each bin, 
and report the equivalent peaking cross-feed yields in Table |l[ Figure |3| shows the data yields after 
subtraction of the peaking BB and cross-feed backgrounds (signal yields). 

The inclusive branching fraction for 6 —s- 57 in each Mfjad bin can be obtained from the partial 
branching fractions by taking into account the fraction, /, of final states that are considered in 
the analysis. The fs are taken from the generator level JETSET fragmentation of the inclusive 
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Figure 3: Signal yields (left plot) and weighted partial branching fractions (right plot), calculated 
with corrected efficiencies. Statistical errors only are shown. 



b ^ s'y signal. Since there are discrepancies in the modeling of our twelve reconstructed final states 
as already noted above, we know that the values of / are not completely reliable. There are also 
uncertainties coming from other final states, e.g. with r] and (p mesons, which have not been looked 
for and therefore cannot be compared with data. 

We divide up the missing final states from our generator-level Monte Carlo into categories 
with a Kl (25%), more than one (4-18%), more than 4-body hadronic systems (0-14%), an 
r]/r]' (5-11%), a sss system (0-7%), and "other" (0-4%), where the missing fractions increase from 
MHad = 1.0 to Mnad = 2.4 GeV/c^. The Kl fraction is known to be 25% from isospin and the 
measured Ks/K~^ ratio in data. For all the other fractions we allow a variation between 0.5 and 1.5 
times the generator fraction, and add these variations in quadrature to give an overall systematic 
uncertainty on the missing fractions. Table ^ summarizes the partial branching fractions with their 
statistical errors, the fraction of studied final states with the corresponding errors, the resulting 
branching fractions with their statistical and systematic errors and the cumulative b ^ sj branching 
fraction, with statistical and systematic errors, in bins of Mnad- 



MHad [GeV/c^ ] 


AB 


•10 


-6 


/(%) 




B 


•10" 






Cumulative 


B -10 


-6 


0.6-0.8 


2 


± 


1 




75 




2.3 


± 


1.3 


± 


0.5 


2.3 


± 


1.3 


± 


0.5 


0.8-1.0 


27 




3 




75 




36 


± 


4 


± 


5 


39 


± 


4 


± 


5 


1.0-1.2 


9 




3 


66 


± 


3 


14 


± 


4 




4 


53 


± 


6 


± 


7 


1.2-1.4 


30 




7 


63 




3 


47 


± 


10 




9 


100 


± 


12 




14 


1.4-1.6 


38 




7 


55 




4 


70 


± 


13 




13 


170 


± 


18 




25 


1.6-1.8 


29 




10 


45 




7 


64 


± 


22 




15 


233 


± 


28 




36 


1.8-2.0 


34 




13 


35 




8 


98 


± 


38 




30 


332 


± 


48 




61 


2.0-2.2 


26 




24 


27 




10 


95 


± 


88 




43 


426 


± 


100 




98 


2.2-2.4 


7 




14 


22 




11 


33 


± 


63 




62 


460 


± 


118 


± 


134 



Table 3: The partial branching fractions (AjB), the fraction of final states that are considered in 
the analysis (/), the total 6 — > 57 branching fractions {B) and the cumulative total b ^ B in 
bins of Mffad- Statistical errors are shown for the A;Ss, and systematic uncertainties for the fs. 
Both statistical and overall systematic errors are shown for the Bs. 
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The systematic errors will be described in Section ^. It can be seen that the branching fraction 
for Mffad < 1-0 GeV/c^ agrees with the exclusive K*^ branching fraction of about 4 • 10^^ and 
that the integral for Muad > 2.0 GeV/c^ is higher than previous measurements of the inclusive 
b ^ sj branching fraction [|l3t 0], although with large uncertainties. Note that these results are 
essentially independent of the modeling of the 6 — > 57 spectrum. 

Figure § shows the branching fractions as a function of Mnad- The same results are shown as 
a function of in Figure ||. 




Figure 4: Branching fraction as a function of Mjjad- The errors are purely statistical. 
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Figure 5: Branching fraction as a function of E^. The errors are purely statistical. 



7 Systematic Errors 

The first category of systematic errors {general systematic) are those which do not depend on the 
hadronic final state. The systematic error on the detection of the high-energy photon is 3.7%, 
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and the systematic error on the total number of B mesons in our data sample is 1.6%. There is 
also a systematic uncertainty of 1% coming from the use of the cos 6t* and Fisher Discriminant 
requirements to suppress continuum backgrounds. These errors give a total contribution of 4.2% 
independent of the hadronic mass bin. 

A second category of errors are due to the fitting procedure and signal definition (fit systematic) . 
These errors are evaluated by varying the fixed parameters used in the fits, and by varying the 
fit procedures. The width of the signal Crystal Ball function is varied by ±0.2 MeV/c^ for both 
signal and peaking backgrounds. This introduces a systematic uncertainty of 5% on the signal yield. 
Varying the Crystal Ball shape to reflect our understanding of the radiative tail introduces a further 
systematic uncertainty of 5%, and varying the peak position introduces a systematic uncertainty 
of 3%. The ARGUS shape parameters are varied within the ranges allowed by the Monte Carlo 
samples. They give a contribution to the systematic uncertainty of 10%, which is limited by the 
Monte Carlo statistics. We also include a systematic error of 20% on the peaking BB yield in the 
region MHad= 1.8-2.4 GeV/c^ due to the method of fitting the BB yield, and to uncertainties in 
the content of the generic BB Monte Carlo simulation. This translates into a 10% systematic error 
in the last three M^ad bins. 

Next we consider errors in the detection efficiencies for the hadronic final states {efficiency 
systematic). We assign systematic errors based on the uncertainties in the corrections to the 
efficiencies that we have applied. The uncertainties coming from the single particle detection 
differences between data and Monte Carlo are between 4.0 and 4.6% as a function of Muad-, and 
the fragmentation differences give uncertainties between 2.1 and 6.3%. 

A final category of errors is due to the generator and modeling we use [generator systematic). 
The errors due to the missing final states not considered in the analysis were discussed in the 
previous section and shown in Table ^. Also, we consider contributions due to changing the modeling 
of the bin 1.0-1.2 GeV/c^ to K*^ and varying the parameters of the inclusive model in our signal 
Monte Carlo. 

The contributions to the systematic errors coming from all sources are summarized in Table 



Muad 


general 


fit 


efficiency 


generator 


total 


[GeV/c2 ] 


systematic (%) 


systematic (%) 


systematic (%) 


systematic (%) 


systematic (%) 


0.6-0.8 


4.2 


17 


4.0 




18 


0.8-1.0 


4.2 


6 


4.0 




8 


1.0-1.2 


4.2 


24 


5.1 


15 


29 


1.2-1.4 


4.2 


12 


7.7 


11 


18 


1.4-1.6 


4.2 


13 


5.5 


12 


19 


1.6-1.8 


4.2 


12 


7.2 


19 


24 


1.8-2.0 


4.2 


17 


5.1 


25 


31 


2.0-2.2 


4.2 


23 


6.4 


38 


45 


2.2-2.4 


4.2 


177 


5.1 


51 


184 



Table 4: All contributions to the systematic error are shown together with the total sum in quadra- 
ture as a function of the hadronic mass interval. 
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8 Fits to the Spectrum 



So far our results have had no significant dependence on the 6 — > 57 model. We now fit the hadronic 
mass spectrum with the shape predicted by the Kagan and Neubert model @. This extrapolates 
from the measured range of Mfjad to give the inclusive branching fraction for h 57. 

In the fits we vary the parameters m;, and Ai, the transition point between the K* and non- 
resonant contributions, and the total normalization, which corresponds to our integrated branching 
fraction. In the plane Ai and ruf,, Figure ^ shows the value of the x^; having fit the transition 
point (resulting around 1.1 GeV/c^) and the normalization for each point. There is a strong 
correlation between these parameters, which is expected in this model, and a shallow minimum at 
rrib = 4.65 GeV/c^ and Ai = —0.48 [GeV/c^]^. Figure 0(a) shows the correlation curve between Ai 
and mi, along which we have the minimum value of ^^'^ Figure 0(b) shows the values of ^is 
a function of irib- Figure 0(c) shows the inclusive branching fraction and the branching fraction for 
Mnad < 2.4 GeV/c^ as a function of rrih, where the result for the minimum value of nib found in 
the free fit is marked with a line. Since our spectrum contains a K* peak with a known branching 
fraction of about 4 • 10~^, a reduction in rub leads to an increase in the non-resonant component of 
the spectrum, and hence to a higher inclusive branching fraction. This increase can be seen both 
in the branching fraction restricted to the measured range, and in the extrapolation over all MHad- 
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Figure 6: Minimum value in the plane of \i,mb. Values of the larger than 5 are set to 5. 
Note the very strong correlation between Ai and mb- 



The free fit to the model parameters only allows us to set an upper bound on mb, since the 
rises rather slowly on the other side of the minimum. We obtain a lower bound on mb by 
measuring the first moment of the distribution, < >, and then extracting A(a^,l/Mg) 
using HQET calculations [p^ p!6[| with the coefficients relevant for our minimum value of E^ = 
2.094 GeV [|l^]. Then, we truncate the HQET expressions to order 0{as,l/M^) and recompute 
mb = Mb — ■^-{cis-, ^° obtain the parameter required in the Kagan and Neubert model. There 

are significant theoretical uncertainties associated with this method coming from varying the scale 
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Figure 7: Results from a free fit to the M^ad spectrum. Plot (a) shows the correlation curve 
between Ai and rrif, along which the is minimized. Plot (b) shows the minimum as a function 
of mi,. Plot (c) shows the inclusive B and the B for Mjjad < 2.4 GeV/c^ as a function of mfo. 
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at which is determined between and 2mfc, from the truncation of higher-order terms, and 

from the dependence of A on the choice of the minimum |16]. 
The value of the first moment, < E^ > |£;^>2.094 GeVi is: 

<E^> |e^>2.094 GeV = 2.35 ± 0.04 (stat) ± 0.04 (syst) GeV, 

which gives: 

A = 0.37 ± 0.09 (stat) ± 0.07 (syst) ± 0.10 (model) GeV/c^. 

where the theoretical uncertainties have been included (model). With the truncation to first order 
m a,, A = 0.49 GeV/c^, which gives: 

mb = 4.79 ± 0.08 (stat) ± 0.10 (syst) ± 0.10 (model) GeV/c^. 



9 Final Results 

We perform final fits to the hadronic mass spectrum with = 4.79 it 0.16 GeV/ taken from the 
first moment analysis. These fits give the inclusive branching fraction: 

B(b S7) = 4.3 ± 0.5 (stat) ± 0.8 (syst) ± 1.3 (model) ■ lO"'^ 

and a range for the parameter Ai: 

Ai = -0.24l°:°3 (stat) ±0.02 (syst)to:ll (model) [GeV/c^]^ 

Figure § shows the superposition of the theoretical spectrum on the data distribution for the 
best fits with mt = 4.79 GeV/c^ and mt = 4.65 GeV/c^. 
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Figure 8: Superposition of the predicted spectrum for mb 
on the observed hadronic mass spectrum. 
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Our values for A, at 0(a1, 1/M^), and Ai, at 0(as, 1/M^), are consistent with other results [18|, 
but our inclusive branching fraction is somewhat higher |p!3|, |l4|. Our error bars are still larger than 
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the previous measurements, but they can be reduced significantly with the addition of more data. 
It is possible to reduce the systematic errors by improvements in the treatment of backgrounds, 
and by studying a larger number of final states to reduce the missing fraction errors. The model 
dependence will also be reduced by fitting the spectrum to determine the model parameters. 
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